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Abstract 

To describe multiple interacting fragmentation continua, we develop a method in which the vibra- 
tional channel functions obey outgoing wave Siegert boundary conditions. This paper demonstrates 
the utility of the Siegert approach, which uses channel energy eigenvalues that possess a negative 
imaginary part. The electron scattering energy in each such channel is rotated upward, giving it an 
equal and opposite imaginary part. This permits a natural inclusion of vibrational continua without 
requiring them to appear as explicit channels in the scattering matrix. Calculations illustrate the 
application of this theory to photoionization, photodissociation, and dissociative recombination. 

PACS numbers: 3.65.Nk,33.80.Eh,34.80.Ht 
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In this Letter, we propose a general method for describing coupling between electronic 
and dissociative continua, based on a Siegert pseudostate basis representation of the vibra- 
tional degree of freedom. The underlying rationale for this idea traces back to the original 
recognition by Kapur and Peierls (1| that the narrow resonances of a scattering spectrum 
can be described in terms of a complex energy eigenstate, with the imaginary part of the en- 
ergy defining a resonance width parameter. This proposal was further developed by Siegert 
in two important ways. First, Siegert demonstrated that a set of eigenstates of any 
Hamiltonian could be chosen to satisfy pure incoming or outgoing wave boundary condi- 
tions in the asymptotic limit. Second, and of importance to our work, Siegert's derivation 
allowed for overlapping resonances of arbitrary width, giving a smooth background term in 
the cross-section. Siegert eigenstates formally correspond to S-matrix poles in the complex 
plane. Sharply resonant features associated with bound states can be identified with poles 
lying on the real axis, while broad background scattering can be described by closely spaced 
eigenstates with finite imaginary parts that serve as a discretized approximation to the true 
continuum. 

Traditionally, the use of Siegert states has been complicated by the nonlinearity of the 
associated eigenvalue problem. Because the wavenumber appears linearly in the boundary 
condition but quadratically (as the energy) in the eigenvalue, the eigenproblem is quadratic, 
and in the past could only be solved iteratively. Tolstikhin et al. |3], [|] recently demonstrated 
how this difficulty may be circumvented using finite range Siegert pseudostates. 

Initially, the true asymptotic boundary condition is replaced by a finite range approxi- 
mation, 



where Rq is a value beyond which the value of the potential is negligible. We seek a solution 
expanded in terms of some primitive basis set 



Here N is the dimension of our basis, and we have selected a non-orthogonal B-spline basis for 
the Uj(R). Inserting this into the Schrodinger equation premultiplying by yf, and employing 
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the boundary value ([!]), we find a matrix equation for the coefficients 

v / 1 rRo 



+ J ° yjl {R)fi[V{r) - E}y 3 (R)dR^j Cj = 0. 
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Note that we have used a Green's theorem identity before substituting in the boundary 
condition, and that the Hamiltonian has been multiplied through by the reduced mass \i. 
Written more concisely in matrix notation, we have a system of the form 

(H-iA;L-£; 2 O)c = 0, (4) 

where Ljji is the surface matrix yj(R )yj>(Ro), Hjji is the matrix 2fxHjj' +yj(Ro)-^-yf(R ), 
and O is the overlap matrix for the spline basis set. 

This equation is manifestly nonlinear, but the method of Tolstikhin et al. allows it to be 
"linearized" by recasting it as a new eigensystem in a basis of doubled dimension. (|| |6|] 
discuss related techniques for solving differential equations where the eigenvalue appears in 
a boundary condition.) We define di = ikci, yielding a trivial second equation ikOc = Od. 
Substituting this into the original eigenequation now gives a linear equation in the doubled 
basis space 

o -o) (cf) lk (-o o ) idj ' ^ 

This is an equation for the eigenvalue A = ik, giving 2N solutions lying either on the (Re 
A)-axis or in conjugate pairs in the right half of the complex plane. 

In their work, Tolstikhin et al. used completeness properties of the Siegert state set to 
construct a Green's function of the Hamiltonian, a scattered solution, and the associated 
scattering matrix, for a variety of single channel model problems. For systems with multiple 
channels, we instead appeal to the well-understood machinery of multichannel quantum 
defect theory (MQDT). 

For resonance series corresponding to high electronically excited intermediates (Rydberg 
states) of diatomic molecules, the most natural description of the system is one with quantum 
defect parameters defined in terms of a fixed internuclear distance R and a well-defined 
projection of the orbital angular momentum A onto the axis of symmetry. This is because 
the electron spends most of its time far from the nuclear core, and when it does penetrate into 



the core, it gains enough speed from falling through the Coulomb potential that the nuclei are 
essentially frozen on the time scale of its motion. The quantum defect functions /xa(-R) in this 
representation, the so-called "body-frame", may either be calculated from highly accurate 
ab initio techniques, or extracted from a semi-empirical fitting of experimental data fl7j. 
In order to connect them with the true asymptotic ionization channels defined in terms of 
Siegert pseudostates of the residual core, j = {v + , iV + }, a frame transformation must be 
performed || [§, where iV + is the ionic rotational momentum, and v + is the vibrational 
quantum number of the pseudostates. In our procedure, we directly evaluate the S-matrix 
by the frame transformation integral 

S jtj , =V(jV + |A) / <P j (R)e 2i7T ^ R ^ j iR)dR(A\N + ') 

The surface term in (6) also arises in the orthonormality relation ||]. A similar transforma- 
tion converts the body- frame transition dipole elements D^(R) into reduced dipole matrix 
elements in the same S-matrix representation, 



^ 5 = (2J+1)^(A|J )( J ) (A|iV+> 

A 
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Here 4>o(R) is the initial vibrational wavefunction, and Jo and J are the total angular 
momenta of the initial and final states of the system, respectively. Note that the Siegert 
pseudostates are never conjugated in these expressions, even when they formally belong to 
the dual ("bra") space. In particular, this means that the quantity labeled as D s ^ below is 
calculated by conjugating only e 47rMA< ^ in the definition above, and not the dipole matrix 
elements directly. 

At this stage of the calculation no information about the long-range behavior of the 
channels has yet been included, and since the body-frame quantum defects are nearly energy 
independent, the resulting S-matrix is typically a smooth and fairly weak function of energy. 
The method of channel elimination [l0| , systematically eliminates flux in all electronic 
channels below the energy threshold for electron escape (the "closed-channel subspace") to 
form a "physical" S-matrix S phys , by taking the proper linear combination of short-range 
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solutions that ensures exponential decay at infinity. For a long-range Coulomb potential, 
this procedure gives 

S^ = S 00 -S 0C (S cc -e- 2 ^)- 1 S C0 . (8) 

Here, (5 is a diagonal matrix of the usual Coulomb long-range phase parameter 7r(i/j) where 
Uj is the (possibly complex) effective quantum number in the jth channel, S is the scattering 
matrix, and the subscripts indicate partitions of the matrices into closed and open subspaces 
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For a Siegert state basis, this physical scattering matrix is in general not unitary, but 
rather subunitary, reflecting the loss of flux at the boundary R via coupling to the Siegert 
pseudo-continuum states. It may be used to calculate partial cross-sections by means of 
conventional formulae, but with the departure from unitarity, 1 — J?. |S^*| 2 , identified as 
the probability |S^T| 2 for scattering into the dissociative continuum. This method also 
provides all quantities necessary to find the partial photoionization cross-section into any 



open channel, af, see Eq. 2.59 of |TT| for further details. The contributions from all open 
channels may then be summed to give the total cross-section for photoionization. 

Alternatively, the total photoabsorption cross-section may be found directly from a "pre- 
convolution" formula first derived by Robicheaux to handle the energy smoothing of densely 
spaced resonances |l2|, |i~3| , 



o-totai (E )= ReD s t 

3(2J + 1) (9) 

[1-Se 2 ^]" 1 [l + Se 2i/3 ] 3 s 

where Re signifies taking the real part of everything that follows, and the f here conjugates 
only the operator, not the entire matrix element. The diagonal matrix written as e~ 2l/3 has 
a nontrivial definition in terms of the quantum defect parameters, it may be approximated 
quite well by taking (3j = ioo for "closed" channels with (E) < Re Ej, and = irvj for "open" 
channels with (E) > ReEj. Here E is the total energy of the system, Ej is the threshold 



energy for channel j, and Uj = l/W2(Ej — E) on the branch where Im v > 0. The utility 
of this expression lies in recognizing that the value of the cross-section at a complex energy 
in the above formula is equivalent to the cross-section at a real energy, smoothed over a 
channel-dependent width Fj = 2Imej. Within the Siegert state formulation, the electron 
energy €j = E — Ej will naturally take on a complex value in any channel where the channel 
eigenenergy Ej is itself complex, while E remains real. 



Given S and e J? either of the two cross-section formulae above may be evaluated, with 
appropriate allowances for the possibility of complex energy eigenvalues. Note that the first 
procedure simply gives a sum over the flux into specific ionization channels, while the second 
gives a single value for the total photoabsorption cross-section. This means that the latter 
will contain information about the solution wavefunction along the R = Ro boundary not 
contained in any of the open ionization channels. In general, the value of atotai will be equal 
to or greater than the sum over the individual <x,-, and any difference may be attributed to 
the effect of coupling to high-lying Siegert states in the continuum. Thus, the difference 
between these two formulae at any energy provides the dissociative cross-section. 

In order to test the validity of this hypothesis, we began by defining a set of Siegert 
pseudostates for the H% internuclear potential. The eigensolutions fall into three classes, 
as shown in Fig. 1. Those lying along the positive (Im A;)-axis are associated with negative 
eigenenergies on the physical sheet of the .E-plane, the bound states of the potential. These 
are the channel thresholds to which the Rydberg autoionization series of the ionization 
spectrum converge, and so we include all of their states. The solutions along the negative 
(Im fc)-axis lie on the unphysical energy sheet, and we reject them as antibound states arising 
from the doubling of the dimension space. The remainder of the solutions fall above and 
below the (Re _E)-axis, corresponding to conjugate solution pairs of the eigenvalue parameter 
A = ik. We select only those with negative Im Ej, ensuring that they obey outgoing wave 
boundary conditions. (This amounts to selecting only those states that contribute to the 
retarded, rather than the advanced, Green's function. See for related discussion.) For 
MQDT matrix elements it is also acceptable to reject states lying very high in the continuum, 
since their Franck-Condon overlap with the bound states is negligible. 

Tolstikhin et al. discuss the unusual completeness relation obeyed by the full set of Siegert 
pseudostates, with an additional factor of 2. Our restricted subset of Siegert pseudostates 
does not, of course, obey that doubled completeness relation. We have confirmed through 
numerical tests, however, that our restricted subset behaves like a complete set, to at least 
10~ 12 accuracy, for representing either L 2 functions confined within the boundary or func- 
tions with purely outgoing wave character at the boundary. For an impressive demonstra- 
tion (in a somewhat different context) of the convergence properties of a similarly truncated 
Siegert basis also used to describe smooth continuum physics, see [ |T5|| . 

In the region of the ungerade H 2 spectrum between 127200 and 127800 cm -1 there are 
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several strongly predissociated resonances, members of the npir,v + = 8 and npir,v + = 
5 series. In each case, our calculated spectrum correctly reproduces them in the total 
absorption cross-section, but shows them as weak or absent in the ionization. Comparisons 
of our results with other theoretical and experimental values [[/], [L6|] for the relative yields 
of selected resonances appear in Table I. Note particularly that we are able to correctly 
describe the strong rotational dependence of the 4pn, v + = 5 branching ratio, a nontrivial 
consequence of subtle channel interactions. 

As a test of the method in an entirely different energy regime we considered the prob- 
lem of dissociative photoionization, a three-body breakup channel accessible only at much 
higher energies. Experimental measures of the ratio between pure ionization and dissocia- 



tive ionization have been performed since the 1970s by a number of researchers [|17], |1| [19 
along with at least one early theoretical calculation [^(J. Since our ionization spectrum is 
a sum over individual channels, we can easily distinguish between contributions from chan- 
nels above and below the dissociative threshold. Our results, plotted against those of past 
experiment and theory, are presented in Fig. 2. 

Finally, we have performed a model calculation demonstrating the utility of our method 
for treating dissociative recombination, particularly in systems where indirect channels 
(those involving scattering into intermediate autodissociating Rydberg states) play an im- 
portant role. Fig. 3 shows the dissociative recombination spectrum of a simplified H 2 model 
potential (neglecting rotation and with R-independent quantum defects), compared with the 
familiar approximation of O'Malley for smooth background scattering by direct processes 



21J| . Our spectrum accurately reproduces this background, and also describes complex in- 



terference effects from the series of resonances converging to each Rydberg threshold. 

Some aspects of the Siegert MQDT method remain poorly understood, and would benefit 
from greater clarification. For example, the utility of a subset of the Siegert basis for MQDT 
depends on ability of that subset to represent all energetically accessible regions of the 
continuum. While this requirement appears from our numerical tests to be reasonably easy 
to satisfy, we have not yet rigorously derived it from the relevant completeness relations. 
Also, it is not presently clear how to extend the energy-smoothed formula to include non- 
Coulombic long-range electronic potentials. 

Other avenues of investigation could provide insight concerning the applicability of our 
method to more complex systems. Polyatomic molecules, for example, might be handled ei- 
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ther by reduction to hyperspherical coordinates [2|[|, or by a multidimensional generalization 
of the Sierget state boundary conditions on an arbitrary hypersurface. Since our method 
yields a full solution to the Schrodinger equation along the boundary at Ro, it should also 
be possible to project onto the continuum functions of different dissociative channels, and 
explicitly resolve partial dissociation cross-sections |p3|1 . Even in its current form, however, 
we believe the Siegert MQDT method offers a simple description of the flux escaping into 
dissociative channels, by working within a channel basis that obeys a physically motivated 
boundary condition. 

This work is supported in part by a grant from the National Science Foundation. We 
thank B. Esry for assistance in the early stages. Discussion with R. Santra, M. Baertschy, 
M. S. Child, C. W. McCurdy, J. M. Hutson, T. N. Rescigno, and B. I. Schneider has also 
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TABLE I. Photoionization and photodissocia- 

tion yields for select ungerade resonances in H 2 
for which the relative yields have been experi- 
mentally observed |16|| . 



State 


Source 


Energy 


% Ion. 


% Diss. 


3pn,v + = 8,R(0) 


Observed 


127248.2 


10(5) 


95(5) 




Theory II 


127246.9 


1 


99 




Present 


127242.2 


1 


99 


5pa, v+ = 4, R(0) 


Observed 


127599.4 


90(10) 


10(10) 




Theory 11 


127602.2 


88 


12 




Present 


127606.8 


76 


24 


4 p7 r, v+ = 5, R(0) 


Observed 


127667.6 


82(5) 


18(5) 




Theory 11 


127665.4 


93 


7 




Present 


127666.6 


97 


3 




Observed 


127599.4 


30(10) 


70(10) 




Theory 


127758.4 


17 


83 




Present 


127759.5 


29 


71 
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Figure 1: Distribution of HjJ" vibrational Siegert pseudostate energies in the complex energy plane 
for angular momentum N + = l. Only the circled states are included as channels in the scattering 
matrix. 




Figure 2: Dissociative photoionization cross-section, as a ratio to the total photoionization cross- 
section. 
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Figure 3: Dissociative recombination cross-section for the model potential, unconvolved (solid) and 
convolved with a Lorentzian of width 0.1 eV (dashed), compared to that resulting from the O'Malley 
formula (dotted). 
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